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Abstract. We model the spin pulse intensity and hard- 
ness ratio profiles of the isolated neutron star RXJ0720.4- 
3125 using XMM-Newton data. The observed variation is 
approximately sinusoidal with a peak-to-peak amplitude 
of 15%, and the hardness ratio is softest slightly before 
flux maximum. By using polar cap models we are able 
to derive maximum polar cap sizes and acceptable view- 
ing geometries. The inferred sizes of the caps turn out to 
be more compatible with a scenario in which the neutron 
star is heated by accretion, and place limits on the mag- 
netic field strength. The hardness ratio modulation can 
then be explained in terms of energy-dependent beaming 
effects, and this constrains the acceptable models of the 
emerging radiation to cases in which softer photons are 
more strongly beamed than harder photons. An alterna- 
tive explanation in terms of spatially variable absorption 
co-rotating in the magnetosphere is also discussed. 

Key words: X-rays: stars, stars: neutron, stars: magnetic 
fields, individual: RXJ0720.4-3125. 



1. Introduction 

From its high ratio of X-ray to optical flux, soft X-ray 
spectrum, 8.4 s X-ray period and location in the Galac- 
tic plane, Haberl et al. (1997) proposed that the ROSAT 
source RXJ0720. 4-3125 is an isolated neutron star (NS). 
A probable optical counterpart has been identified (Motch 
& Ha berl |1998j ; Kulkarni & van Kerkwijk |1998j ). Paerels 
et al. Q2001| , Paper 1) present spectra of RXJ0720.4-3125 
using XMM-Newton, and find that there is no evidence for 
absorption lines and edges in the X-ray spectrum. They 
are able to refine the ROSAT spectral fits using black- 
body spectra of 86.2±0.3 eV absorbed by a cold absorber 
of 6.0 ±0.4 x 10 19 cm" 2 . 

Paper 1 mainly addressed the spectral characteristics 
and the behaviour of RXJ0720. 4-3125 averaged over the 
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8.4 s spin period. In this paper we examine the implica- 
tions of the phase-resolved data. 

2. Observations 

Observations of RXJ07 20. 4- 3125 were made using XMM- 
Newton (Jansen et al. |200l| ) on 2000 May 13 (orbit 78). 
For the phase resolved data we have concentrated on the 
EPIGPN data, as the EPIGMOS provides insufficient 
time resolution in the modes selected for the observa- 
tion. The raw data were processed using the version of 
the XM M-Newton Science Analysis System (Watson et 
al. [2001D released on 2000 July 12. 

The particle background was significantly higher at the 
end of the observation, so those data were removed from 
the analysis, leaving a total observing time of 49 ks. Light 
curves were extracted from the EPIC-PN camera using 
an aperture ~ 30 in radius centered on RX 0720.4-3125, 
chosen so that the aperture did not cover more than one 
CCD. This radius encompasses ~ 90% of the integrated 
PSF (Aschenbach et al. 200d| ). We accumulated counts 
over the energy ranges 0.1-0.4, 0.4-0.8 and 0.8-1.2 keV 
and phase-averaged them over 42 phase intervals on the 
8.391 s period of Haberl et al. fll997| ). The spin profile in 
the band 0.1 — 1.2 keV is shown in Figure |[ together with 
the hardness ratio (soft/medium) variation and the best 
fit discussed in Section |4j 

2.1. Phase-folded X-ray light curves 

The general shape of the phased intensity curve is similar 
to that in Haberl et al. ( |1997| ) and is described in Paper 1. 
It appears approximately symmetrical and sinusoidal with 
an amplitude ~ 15%. The hardness ratio is also seen to 
vary: it is softest around flux maximum and the amplitude 
is smaller at ~ 10%. The phasing of the hardness ratio 
curve is slightly but significantly earlier than the intensity 
curve. A cross-correlation of the two curves indicates that 
the phase difference is —0.048. 

The most likely source of the flux variation is the 
changing visi bility of the heated magnetic polar caps 
(Haberl et al. 1997 ) with the rotation of the NS. Assum- 
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Fig. 1. The total flux (top) from EPIC- PN ph ase aver- 
aged on the 8.391 s period of Haberl et al. ( 1997[ ) with the 
(0.1 to 0.4)/(0.4 to 0.8) keV hardness ratio (bottom). The 
flux is normalised to the model flux at phase 0.0. Also 
shown is the best fit to the data solid polar cap model 
with beaming and temperature variation (see Section ^ 
for details). 



ing symmetry, it is possible that the rotation period is 
twice that deduced by Haberl et al. ( 1997 ), depending on 
whether one cap or two are seen each rotation. A Fourier 
analysis excludes any harmonic or subharmonic to the 
level of ~ 10% of the 8.391 s period but it is possible 
that odd-even effects may be present at lower levels. We 
have therefore also folded the data on a 16.782 s period. 
This gives rise to a lightcurve with two peaks per spin pe- 
riod with a small difference in amplitude. In this case the 
hardness ratio is also found to show two peaks per spin 
period. 



2.2. Phase resolved spectra 

We have selected data between phases 0.85-0.05 and 0.15- 
0.35 to examine more closely the change in spectral param- 
eters in the "soft" and "hard" phase, i.e. at the maximum 
and minimum of the lightcurve, respectively. Spectra were 
background corrected using the same CCD on which the 
source was detected. In our fits we used the response file 
epn_fs20_sY9_thin.rmf produced by one of us (FH, version 
2000 July) and we selected only single pixel events. 

As in Paper 1 we have used a simple absorbed black 
body model in XSPEC. The spectral fit suggests that the 
majority of the variation is caused by a change in ab- 
sorption, rather than a change in temperature (Figure ^). 
The fits are consistent with a temperature of 85 eV, with 
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Fig. 2. Confidence contours in the temperature-absorbing 
column plane for the phases at the extremes of the hard- 
ness ratio variation in Figure [l]. Contours are for 68%, 90% 
and 99% confidence level. 



the measured temperature difference between flux maxi- 
mum and minimum of less than leV. On the other hand, 
the best fit of the "hard" phase requires a doubling in 
absorption from 4 x 10 19 cm~ 2 to 8 x 10 19 cm~ 2 . These 
parameters are consistent with the phase-averaged values 
from RGS data derived in Paper 1. While residual uncer- 
tainties in the EPIC-PN calibration at these soft energies 
may introduce some uncertainties into the absolute values 
of the derived absorptions (see Paper 1), it is unlikely that 
the relative change can be discounted. 



3. Analysis of the light curve 

3.1. Polar cap models 

The high S/N ratio of the X-ray pulsation in the XMM- 
Newton data provides the opportunity for a more de- 
tailed modelling than was possible using the ROSAT data. 
In this section, we first present the best fit of the light 
curve only, re-adapting the formalism originally derived 
by Pechenick et al. (1983). Assuming that the intensity 
variation is due to the changing visibility of heated circu- 
lar magnetic polecaps and including gravity effects, those 
authors were able to derive a semi-analytical expression 
for the observed brightness, in terms of elliptic integrals. 
Their final result depends on four variables: the angular 
radius of the polecap a, the angle between the dipole and 
rotation axes (3, the angle between the line of sight and 
the rotation axis 7, and R/2M, where M and R are the 
mass and radius of the NS (G = c = 1). In addition, a 
prescription for the beaming function, /(£), and a phase 
angle are required to complete the viewing angle calcula- 
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f(S) 


R/2M 




1 


4 


53° 




3 


44° 




2.5 


26° 


cos 8 


4 


63° 




3 


61° 




2.5 


59° 


exp(— 5) 


4 


65° 




3 


63° 




2.5 


61° 



Table 1. The maximum angular size of the polar caps for 
different beaming prescriptions and compactness parame- 
ters. 

tion. The equations are symmetric for (3 and 7 and R/2M 
ranges between 2.5 and 4.0 for standard NS masses and 
equations of state. 

Using this formalism, we have calculated model light 
curves for a range of parameter values. We have then used 
a Levenberg-Marquardt algorithm to fit the phased data 
by least squares. In order to maximise the precision of the 
fit we have used the spectral range of 0.1 to 1.2 keV and 
we have weighted the fit using the reciprocal of the squares 
of the uncertainties. 




Fig. 3. The xt over the (/3, 7) plane for a cap angular ra- 
dius a = 10° and R/2M = 3.0. Here f{6) = 1. Contours 
are for ^=1.0 and then the 68% and 90% confidence in- 
tervals. The arc-shaped region and the small semicircular 
regions near the axes correspond to periods of 8.391 s and 
16.782 s, respectively. 



3.2. Fit to the models 

In the case of isotropic emission (no beaming, f(S) = 1) 
the best fit to the data occurs for a polar cap radius of a = 
35°, dipole and viewing angles of (3 = 9.3° and 7 = 42.7°, 
and R/2M = 4.2. The quality of the fit is good, with 
X 2 =0.96, and for these values the cap at the other pole is 
at most barely visible. The values of a and (5 are the most 
sensitive to the angular distribution of the radiation field: 
when we change the beaming prescription by adopting 
a "pencil" model (f(S) = cos 5), the best fit occurs for 
a = 12.4, $ = 4.2, 7 = 41.3 and R/M = 4.2 ( X l=l.U). 

However, the quasi-sinusoidal nature of the phased in- 
tensity data has the result that the phase space is de- 
generate for a wide range of input parameters. We have 
therefore taken the approach of fixing the cap radius and 
R/2M and exploring the (/3, 7) plane for good fits. In the 
case of isotropic emission we have found that, for most 
cap radii and for 2.5 < R/2M < 4.0, the acceptable pa- 
rameter space is concentrated in a hyperbola described 
approximately by f3 ~ C/j, where C ~ 400 for angles 
in degrees (see figure ||). As the cap size increases the 
best fit moves from the center of the arc where ft = 7 to 
(/3, 7) ~ (10,48), and the constant C increases slightly. 

A physically more important result is that, for a given 
value of R/2M, the fits to the data constrain the max- 
imum angular radius of the cap. By increasing a, the 
xt of the best fit in (/3, 7) shows a sharp increase be- 
yond a threshold value. We can therefore identify the 
maximum polecap radius beyond which a best fit can be 
rejected at the 90% confidence level (x 2 =4.61). Results 



are summarized in table 1 for isotropic emission, a pencil 
model and a very extreme case in which the angular distri- 
bution of radiation decays exponentially with S. Smaller 
R/2M require smaller maximum polecap radii because of 
the increasing bending of the light ray trajectories from 
stars with large compactness parameter. In the presence 
of strong beaming effects, larger caps are allowed, but, at 
least when using simple, non energy-dependent beaming 
prescriptions, a polecap larger than ~ 60° — 65° can be 
rejected at a confidence level of 90%. 

3.3. Constraints on B from the size of the polar caps 

Independently on the interpretation of RXJ0720. 4-3125 
as a cooling or an accreting NS, the existence of an upper 
limit for the cap size requires the NS to be magnetised. 

A cooling NS interpretation implies magnetic fields 
much larger than in the accreting case. Such a high field, 
typically 10 12 — 10 14 G, yields temperature variations 
with the magnetic colatitude, x, and modifies the emer- 
gent spectrum as function of viewing angle (Greenstein 
& Hartke 1983. Page 1995). A simple expression for the 
surface temperature profile is (Possenti et al. 1996): 



T, 



surf 



eff 



[1 - 0.47(1 - K)] 



1/4 



K + (4 - K) cos 2 x 



1 + 3 cos 2 x 



1/4 



(1) 



where K = k±/ku is the ratio of the coefficients of ther- 
mal conductivity orthogonal and parallel to the field: 
K - 10~ 2 - 10~ 4 for B - 10 11 - 10 13 G. Not even the 
maximum values of the polar cap size derived in Section 
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|3.2| for strong beaming are fully compatible with this kind 
of interpretation: the temperature gradient given by equa- 
tion (|l|) is so smooth that, even in the (unrealistic) limit 
case K = the associated blackbody luminosity drops 
only by a factor 1 /2 for a ~ 60° and by an order of mag- 
nitude for a ~ 77°. Both the emergent spectrum and the 
degree of modulation of the signal are not strongly de- 
pendent on K. Despite the real angular distribution of 
the specific intensity being unknown and the fact that a 
proper modelling of the relative brightness should include 
the temperature profile, it seems difficult to explain the 
(relatively) small hot spot inferred from our calculations 
in terms of a surface temperature variation during the NS 
cooling phase. 

If the polar caps are heated by accretion, the incoming 
material is channeled to the poles only by those field lines 
which have passed the Alfven radius. In the case of an 
aligned, dipolar B-field this gives r cap ~ i? 3 / 2 /^ 2 , where 
reap, r A are the cap and the Alfven radius, respectively 
(see e.g. equation (9) in Treves et al. 2000D - Thus, the 
previous values of the maximum polecap radii translate 
for the case of isotropic emission to B > (5xl0 5 ,6x 

10 4 ,2 x 10 4 )Afif G, where M n = A//(10 11 gs- 1 ). 

The minimum polecap radius is not constrained from 
the data: we find that good fits can be obtained even for 
a = 1° for all 2.5 < R/2M < 4.0. In the framework 
of the accreting scenario, the existence of an upper limit 
B ~ 10 10 G can be used to constrain the minimum a: for 
larger fields, in fact, accretion will be completely hindered. 
This gives a minimum a ~ 1°, which is indeed consistent 
with the data. Such low field strengths may provide one 
of the first indication that a significant field decay had oc- 
curred in an isolated NS, or, as it has been also suggested, 
may indicate that this object is the outcome of the com- 
mon envelope evolution of a binary system (see Treves et 
al. 2000 and references therein, Wang 1997). 



4. Simultaneous analysis of light curve and 
hardness ratio variation 



The Pechenick et al. (1983) formalism assumes emission 
from a uniform temperature polecap, and it does not pre- 
dict any hardness ratio changes through the spin cycle. 
However, zones of different temperatures can be used in 
conjunction with the spectral response of EPIC-PN and a 
spectral model for the emission to calculate the observed 
fluxes and hardness ratios. We use a blackbody to simu- 
late the thermal emission as this provides a good approx- 
imation to the flux distribution of more detailed atmo- 
spheric models (at least in the X-ray band, see Zavlin et 
al. 1996, Zane et al. 2000). However the angular distri- 
bution of the emerging radiation in these models show a 
variety of beaming profiles, depending on chemical com- 
position and magnetic field, which must be accounted for 
separately in the modelling, as implemented in Section 3. 



Model 


8.391 s period 


16.782 s period 


2-polc 

a < 35° 


2-polc 

a > 35° 


1-polc 


2-polc 

a < 35° 


2-polc 

a > 35° 


1-polc 


2 polar caps at dif- 
ferent temperatures 


x A 


x F 


x c 


X D 


x B 


x B 


2 polar caps with 
different absorptions 


x A 


x F 


x c 


x D 


x B 


x B 


2 equal polar caps 
with gradient in ab- 
sorption 


x A 


O 


O 


o 


x B 


x B 


2 equal polar caps 
with gradient in 
temperature 


x A 


x B 


x B 




x B 


x B 


2 equal polar caps 
with beaming 


x A 


m 


m 


• 


x B 


x B 



Table 2. Summary of acceptable models: 
• Acceptable fits possible. 

O Probably unacceptable on grounds of poor model fit 
and unphysicality of absorption model 

A Excluded on the grounds of not being consistent with 
the 8.391 s period. 

B Excluded on the grounds of not being consistent with 
the 16.782 s period. 

C Excluded on the grounds that it predicts no hardness 
ratio variation. 

D Excluded on the grounds that the hardness ratio vari- 
ation would be on a 16.782 s period rather than the 
8.391 s observed. 

E Excluded on the grounds that the temperature differ- 
ence required to produce the hardness ratio changes is 
not consistent with the spectral fit. 

F Excluded on the grounds of the shape of the hardness 
ratio variation. 

The two main observational features are the simultane- 
ous variations of hardness ratio and luminosity and their 
shift in phase. In addition, the spectral fit at the minimum 
of the lightcurve implies a larger absorption. By working 
on the premise of minimizing the number of additional 
assumptions, we investigated the parameter space using 
different simple models for the emission. We allow for a 
temperature variation inside the caps, for beaming effects, 
for variations in absorption and for asymmetric conditions 
in the two spots. Despite the temperature of the dominant 
thermal component not varying during a period, the pres- 
ence of a second component with 60 eV < T < 90 eV (or, 
equivalently, of a smooth gradient of T) is still acceptable 
on spectral basis. 

Our findings are summarized in Table 2. With "1-pole" 
and "2-pole" we refer to models in which one or two polar 
caps become visible during the rotation. Quite interest- 
ingly, most models can be discarded since they do not 
allow adequate fits. This is denoted by an x in the ta- 
ble, together with the reason for their exclusion. We find 
that limb darkening variations are the best candidate to 
give rise to the expected modulation. If the emission is 
isotropic, the fact that the hardness ratio is softest near 
the flux maximum can be fitted only if the temperature de- 
creases toward the center of the cap; however, the spectral 
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dependence is inconsistent with the observed spectrum. 
On the other hand, by using different prescriptions for 
f(5), we can also allow for a more realistic temperature 
profile, increasing towards the centre. The observational 
features require a beaming in the softer component much 
stronger than at highest energies. If this is the case, the 
change in Nh found with the simple absorbed blackbody 
model in XSPEC should mimic artificially the soft ex- 
cess/deficit at different phases: this is naturally predicted 
by this kind of angular distribution. Our best fit model, 
X^=l-23, is shown in Figure [l]. The fit to the intensity 
data is good, but is not optimised: because of the number 
of degrees of freedom (now including the choice of spectral 
parameters and phase offsets) a fit was sought manually to 
the intensity data which would be both sufficient in formal 
terms and predict adequately the hardness ratio variation 
in mean value, amplitude and phase. In this model, we 
allow for a cooler ring with T = 70 eV, around a central, 
hotter region of T = 90 eV; the phase offset between them 
is 0.113. The value of the column density has been fixed 
at 6 x 10 19 cm~ 2 , to be consistent with Paper 1. The an- 
gular sizes of the soft and hard components are fitted by 
ot-hot = 33°, a co id — 38°, while the other parameters are 
09, 7, R/2M) = (5.3°, 43.5°, 4.2). The angular dependence 
of the emerging radiation is assumed to change from a pen- 
cil model, f(5) = cos 6, at low temperatures to isotropic 
emission, /(<$) = 1, at high temperatures. No reasonable 
fit can be found when the beaming is stronger at high en- 
ergies, unless we require the polar cap to be colder at the 
center. 

The inferred behaviour of the degree of anisotropy is 
contrary to what is physically expected at the simplest 
level in both cooling and accretion atmospheres. In an 
atmosphere in LTE where the main radiative process is 
bremsstrahlung, the opacity decreases with the photon 
energy. This, in turn, implies the anisotropy in the an- 
gular distribution being smaller at low frequencies (since 
these photons escape from the outermost atmospheric lay- 
ers) and increases with the energy when photons decouple 
in the deepest, h ottest layers. Moreover, as pointed out 
by Zavlin et al. (1996) for pure H cooling models, the 
maximum of the intensity spectrum shifts towards lower 
energies with increasing 5, contrary to the hardness ratio 
variations. In order to match the hardness ratio variation, 
the energy-dependence of the degree of anisotropy must 
be reversed. For hydrogen atmospheres, that may happen 
only in the very high energy part of the spectrum, were the 
angular dependence tends to match the Ambartsumian- 
Chandrasekhar function for scattering dominated propa- 



gation (see Figure 8 in Zavlin et al. 1996). 

The observed change in Njj may be real (see Paper 1), 
but, at least within our simple assumptions, models with 
different levels of absorption do not produce acceptable 
fits. Moreover, a variation in absorption over the spin cy- 
cle require the absorbing material to be close to the NS, 
probably co-rotating in the magnetosphere. In this case 



the light cylinder limits the maximum distance to the ab- 
sorber to Rmax ~ 4x 10 10 cm. On the other hand, in order 
for the material to be substantially neutral, the ionisation 
parameter £ = L/nR 2 must be less than 1. This in turn 
requires a particle density n > 10 10 cm~ 3 , for R < R max 



and a luminosity of 10 erg s (Haberl et al. 1997). A 
similar result can be obtained from the inferred increase 
in absorption of 4 x 10 19 cm -2 . For this to occur within 
Rmax then n > 10 9 cm~ 3 . For comparison, typical den- 
sities of the interstellar medium are n ~ 1 cm -3 . The 
implication is that a reservoir of accreting material is held 
in the magnetosphere, in regions where gravitational, cen- 
trifugal and magnetic forces balance and confine it. These 
densities may be considered uncomfortably high, although 
we are unable to exclude them on current knowledge. We 
return to this in Section 5. 

5. Conclusions 

By exploiting the superior quality of the data obtained 
from XMM-Newton , we have been able to present a de- 
tailed analysis of the spin pulse profile of RXJ0720. 4-3125. 
We have modelled these data by usin g polar cap mod- 
els based on the Pechenick et al. ( 1983 ) formulation. The 
smooth nature of the data does not allow a unique conclu- 
sion about the best fit values to be reached, but neverthe- 
less we are able to exclude a large volume of the parameter 
space and to derive an upper limit for the size of the emit- 
ting region. The derived values are at best marginally com- 
patible with an interpretation in terms of surface temper- 
ature variation during the cooling phase of a middle-age 
NS. An accretor model is more satisfactory in accounting 
for the intensity plus hardness-ratio variations. Accreting 
NS appear to be relatively rare objects, as recently found 
by Popov et al. |2000j fr om a statistical evolutionary anal- 
ysis. RXJ0720. 4-3125 may therefore be more important 
than perhaps initially appreciated. 

We have investigated a number of possible explana- 
tions for the intensity and spectral variations over the 
spin period. The long period of RXJ0720. 4-3125 is some- 
what unusual and hints at a possible evolutionary link 
between dim NS, soft gamma-ray repeaters (SGRs) and 
anomalous X-ray pulsars (AXPs) (Haberl et al. [1997 , Al- 
par 2000 ) , whereby isolated NS in the propellor phase are 



the progenitors of AXPs and SGRs. If the hardness ratio 
variation is indeed caused by absorption, then the large 
amount of material close by to the neutron star may be 
explained in terms of matter being propelled outwards. 
However the luminosity in this phase is caused by en- 
ergy dissipation in the neutron star, and we have already 
found (Section 3.3) this requires polar cap sizes which are 
uncomfortably large. Moreover, such an evolutionary link 
would produce significantly larger numbers of AXPs and 
SGRs (Mereghetti, private communication). 

The most plausible interpretation is that the observed 
modulation is produced by beaming effects in the emerg- 
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ing radiation. The angular distribution of the emerging 
radiation is still a matter for further investigation, and 
is extremely sensitive to details in the radiative computa- 
tions and to a number of assumptions about chemical com- 
position and magnetic field (see e.g. Zavlin et al. 1996| ), 



Rajagopal et al. 1997). Our modelling provides observa- 
tional inputs for these investigations by suggesting that 
the softer X-rays from the heated cap are more strongly 
beamed than the harder X-rays. 
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